Cavity-enhanced light scattering in optical lattices to probe atomic quantum statistics 
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Different quantum states of atoms in optical lattices can be nondestructively monitored by off- 
resonant collective light scattering into a cavity. Angle resolved measurements of photon number and 
variance give information about atom-number fluctuations and pair correlations without single-site 
access. Observation at angles of diffraction minima provides information on quantum fluctuations 
insensitive to classical noise. For transverse probing, no photon is scattered into a cavity from a 
Mott insulator phase, while the photon number is proportional to the atom number for a superfluid. 
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Studies of ultracold atoms in optical lattices link vari- 
ous disciplines. Fundamental quantum many-body theo- 
ries, formulated initially for condensed matter, can be 
tested in better controllable atomic systems e.g., 
strongly correlated phases, quantum simulators. Such 
studies influence different areas [l| : quantum information 
processing, ultracold collisions, exotic molecules, etc. 

While mean-field approaches describe only the average 
atomic density, the main goal is to study quantum prop- 
erties of these gases. They are most prominent in lattices, 
where one has phase transitions between states of similar 
density but radically different quantum fluctuations. 

Standard methods to measure quantum properties are 
based on matter- wave interference of atoms released from 
a trap [3] destroying the system. "Bragg spectroscopy" 
using stimulated matter-wave scattering by laser pulses 
proved successful 0, 0] but destructive. Alternative less 
destructive methods observing scattered light were pro- 
posed mainly for homogeneous Bose-Einstein conden- 
sates (BEC) 0,11, Hill, but not yet implemented. 

Here we show that specifically for periodic lattices, 
light scattering can help to overcome experimental dif- 
ficulties. In contrast to homogeneous BECs, scattering 
from a lattice allows to determine local and nonlocal 
correlations without single-atom optical access using the 
suppression of strong classical scattering at Bragg min- 
ima and monitoring much richer angular distributions. 
This looks extremely useful for studying phase transi- 
tions between, e.g., Mott insulator (MI) and superfluid 
(SF) states, without destruction, since various quantum 
phases show even qualitatively distinct scattering. 

Joining two fields, cavity quantum electrodynamics 
(QED) and ultracold gases, will enable new investiga- 
tions of both light and matter at ultimate quantum levels, 
which only recently became experimentally possible ^] . 

Our model is based on nonresonant interaction, not re- 
lying on a particle level structure. Thus it also applies to 
molecular physics, where new quantum phases were ob- 
tained • It can be also applied for semiconductors 11 1 , 
as, e.g., were used for BEC of exciton-polaritons |12|. 

Model. We consider N two-level atoms in an optical 
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FIG. 1: (Color online) Setup. A lattice is illuminated by a 
probe at the angle So which is scattered into a cavity at 6i. 



lattice with M sites. A region of K < M sites is il- 
luminated by probe light which is scattered into another 
mode (cf. Fig.[T]). Although, each mode could be a freely 
propagating field, we will consider cavity modes whose 
geometries (i.e. axis directions or wavelengths) can be 
varied. A related manybody Hamiltonian is given by 
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where ao (oi) are the annihilation operators of the probe 
(scattered) light with the frequencies uJo,i, wave vectors 
ko,i, and mode functions Mo,i(r); ^'(r) is the atom-field 
operator. In the effective single-atom Hamiltonian Hai, 
p and r are the momentum and position operators of an 
atom of mass rua trapped in a classical potential V^i(r); 
go is the atom-light coupling constant. We consider the 
field-atom detunings Aia = wi —Wa larger than the spon- 
taneous emission rate and Rabi frequencies. Thus, in Hai 
the adiabatic elimination of the upper state was used. 

Assuming weak fields ao,i, we expand ^'(r) in Eq. ^ 
using localized Wannier functions corresponding to Vci(r) 
and keep only the lowest vibrational state at each site: 
^('') = '^iLi^i'^i''^ ~ ''^i): where bi is the atom annihi- 
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lation operator at site with coordinate r^. Substituting 
this in Eq. ([T L one can get a generahzed Bose-Hubbard 
Hamiltonian l| including Hght scattering. However, in 
contrast to our previous work and "Bragg spec- 
troscopy" [1], we do not consider lattice excitations here 
and focus on scattering from atoms in a prescribed state. 
Neglecting atomic tunneling, the Hamiltonian reads: 
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where fii = b\bi. We define the operator of the atom 
number at illuminated sites as Nk = X^^^i ^i- I^^r a deep 



lattice the coefficients JJ 
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J drw'^ir - r,)w,*(r)u„(r) 
reduce to J-™ = u*{ri)um{y^i) neglecting atom spreading, 
which can be studied even by classical scattering [f 4!] . 

The Heisenberg equation for the scattered light in the 
frame rotating with loq (Aqi = luq — lui) thus reads: 
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where k is the cavity decay rate and ao will be assumed 
a classical field given by a c-number constant. 

Light properties. Though the dispersion shift of a 
cavity mode is sensitive to atom statistics through hi , we 
assume it is much smaller than k or Aqi. A stationary 
solution of Eq. ^ for oi and photon number then reads 



ai = CD, 
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with C = i(?Qao/[Aoa(iAoi — k)] and the coefficients 
^4(^*0,^1) = u\{Yi)u{){Yi). This expression of the hght 
operators through the atomic ones is a central result here. 

For a ID lattice of period d and atoms trapped at Xm — 
md (to = 1, 2, M) the mode functions are Mo,i(rm) — 
exp(«TOfco,ia;rf) for traveling and Mo,i(r,„) = cos{'mkQ^ixd) 
for standing waves with ko^x = |ko^i| sin0o,i (cf. Fig.[T]). 
For the atomic quantum state we use the assumptions: 
(i) the mean atom number at all sites is {hi) = n = N/M 
{{Nk) — Nk = riK) and (ii) the pair correlations {hihj) 
are identical for any i ^ j, which is valid for a deep 
lattice, and will be denoted as {hah^) (with a ^ b). 

Thus, (ai) ^ {D) = — "-^ showing that 

the field amplitude only depends on the mean density 
and exhibits the ai^ular distribution of classical diffrac- 
tion A{9o, 61) = J2i=i ^j(^o, ^1) with diffraction maxima 
and minima. The central point now is that the photon 
number ^ is not just the amplitude squared, but we get 
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TABLE I: Statistical quantities of typical atomic states. 
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where Shi = hi — n giving {ShaShf,) = {hah},) — v?, and 
{Sn^) equal to the variance (Ari^)^ = {nf) — -n? . Thus, 
the intensity is sensitive to atomic quantum statistics 
via the density-density correlations {hihj) different for 
particular states. Besides the classical angle dependence 
Ap, the second term in Eq. (|4aP reflects fluctuations and 
has a completely different dependence. Particularly in a 
lattice, scattering is sensitive not only to the periodic 
density, but also to periodic fluctuations, leading to the 
observable difference between states with and without 
nonlocal correlations. Analysis of quadrature variances 
gives results similar to analysis of the noise quantity R. 

For two traveling waves, Eq. (|4a|) gives the structure 
factor considered in Ref. |7j] on homogeneous BECs. We 
show that a more general form including standing waves 
gives new measurable quantities beyond structure factor. 

The intensity fluctuations of the scattered light de- 
pend on the fourth moments of the atomic number op- 
erators and four-point density correlations {hihjhkhi). 
For example, the photon-number variance is given by 
(Anph)^ = {nl^) ~ (nph)2 = \C\H{\D\^) {{0^ + 

\C\H\D\^). 

To discuss examples of different scattering we summa- 
rize statistical properties of typical states in Table HI For 
light scattering, the most classical state corresponding to 
pointlike atoms is MI. Here the atom number at each site 
hi does not fluctuate and we have no pair correlations. 
Hence we see from Eq. (|4a| that the zeros of classical 
diffraction [^(6*0, 6*1) = 0] are zeros of light intensity. 

This is different for a SF where each atom is delocalized 
over all sites leading to number fluctuations at a given 
site and a.t K < M sites; the atoms at different sites are 
anticorrelated. At a classical diffraction zero we still find 
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a photon number proportional to the atom number N. 

A coherent state approximates a SF but without pair 
correlations. In the limit N, M ^ oo, it well describes 
scattering from a small region {K 4; M) of a partially 
illuminated superfluid (SF/f). However, we proved that 
even in this limit it fails to describe scattering at angles 
of Bragg maxima from a large lattice region {K ~ M). 

Example. Let us now show the most striking predic- 
tions of this model at the basic example of a probe trans- 
verse to the lattice {Oq = cf. Fig.[T]). The scattered light 
is collected in a cavity along the lattice {9i = 7r/2) with 
atoms trapped at the antinodes {d — A/2) [3, [3]. 

([3]) here gives al- 
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most zero average field amplitude independently on the 
atomic state. This reflects the opposite phase of light 
scattered from atoms separated by A/2 (diffraction min- 
imum). However, the cavity photon- number is propor- 
tional to {D*D) = ((^2} - {hanb))K [cf. Eq. (Ha])], which 
is determined by statistics of a particular state. Thus, 
atoms in a MI state scatter no photons, while a SF scat- 
ters number of photons proportional to the atom number: 

(ai)Mi (ai)sF = 0, but 
{a\ai)Mi = 0, (a|ai)sF \C\'^Nk. 

Moreover, the photon number fluctuations (Anph)^ 
are also different for various atomic states. In the MI 
state, the variance (A|D|2)2jj = {\D\'^)mi- {\D\'^) 
whereas in SF, there is a strong noise {A\D\'^)gp 
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Coupled light-matter dynamics in a cavity can lead to 
a new self-organized phase 15| with atoms trapped at 
every second site (d = A), which gives D = X]fc=i ~ 
Nk O. If this state is a MI with d — \, the number 
of photons (a|ai)scif-org — \C\'^N'^ is proportional to the 
atom number squared and has a superradiant character. 

Angular distributions. We will quantitatively dis- 
cuss angular intensity distributions for scattering be- 
tween two traveling waves, where Eq. (j4bp reduces to 



R = (ShaShb) 



sin^ {Ka^/2) 
sin^ {a- 12) 



{{Sn") - {Sha5m))K. (5) 



While \A\'^ in the first term reproduces classical diffrac- 
tion with a_ = koxdsinOo — kixdsm9i, the second term 
in Eq. (|4b|l is simply isotropic. Thus, the noise quan- 
tity is zero for MI, Rmi = 0, nonzero but isotropic for 
the coherent state, i?Coh = Nk, and angle dependent for 
a SF. In a SF, even small pair correlations {SfiaShb) = 
—N/M"^ give a large contribution near diffraction max- 
ima (a_ ~ 2ttI, / = 0, 1, ..), where the geometric factor is 
K^, invalidating the coherent-state approximation. 

Figure [2] displays those angular distributions. Classi- 
cal diffraction 1(13)1^ with the only possible zero-order 
maxima at 9i — 0,Tr (d = Ao,i/2, = 0) is shown in 
Figl^Ia). R for the coherent and SFk states are plotted 
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FIG. 2: (Color online) Intensity angular distributions for two 
traveling waves, (a) Intensity of classical diffraction; (b) noise 
quantity R ((5]) for coherent atomic state (constant 1, line A), 
SF with all sites illuminated K — M (curve B), and MI state 
(constant 0, line C); (c) the same as in (b) but for partially 
illuminated SF with K = M/2. iV = M = 30, 0o = 0. 
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FIG. 3: (Color online) Intensity angular distributions for two 
standing-wave modes, (a) Intensity of classical diffraction; 
(b) noise quantity for coherent state; and (c) for SF (curve 
A) and MI (constant 0, line B). N = M ^ K ^ 30, 6*0 = O.Itt. 



in Figs. IDJb), (c). For MI, R = 0. In SF, there is a 
noise suppression at maxima, which is total for all sites 
illuminated, K = M, and partial for K = M/2. 

In a maximum, _D ([3]), is reduced to Nk- Thus, the 
field amplitude is determined by Nk = nK, the intensity 
depends on {D*D) = {N'^), while R = (AiV^)^ gives the 
atom-number variance at K sites, which refiects the total 
and partial noise suppression in Figs.[2fb) and[2jc), since 
(Nk) fiuctuates for K < M. In diffraction minima, the 
field is zero, but the intensity is proportional to (n^) — 
(fiafib)- Under scattering of spatially incoherent light, 
the intensity is isotropic and proportional to (n^). 

So, in optical experiments, varying the geometry, the 
global statistics of K sites, local single-site statistics, and 
pair correlations can be obtained even without a single- 
site access. Thus, light scattering gives a way to distin- 
guish between atomic states. As shown by Eq. ^ and 
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Fig.[2l MI and SFm states are different in diffraction min- 
ima and in incoherent light. They are indistinguishable 
in maxima. SFm and coherent states differ in maxima 
only. The MI and coherent state are different at any 
angles. 

The noise quantity or photon statistics are different in 
orders of Nk for various states. Nevertheless, for large 
Nk, there could be practical problems to subtract large 
values in a maximum. In Refs. 6], this even led to a 
conclusion about the state indistinguishability by inten- 
sity measurements. In contrast to homogeneous BECs, 
in lattices, this problem has a natural solution: measure- 
ments outside maxima are free of the strong classical-like 
part of scattering and thus directly reflect fluctuations. 

A classical analogy of different light scattering consists 
in different density fluctuations. A quantum treatment 
gives a deeper insight. Superfluid state is a superposition 
of all possible multisite Fock states giving distributions 
of N atoms at M sites. Various Fock states become en- 
tangled to scattered light of different phases and ampli- 
tudes. In contrast to a classical case (and MI with the 
only multisite Fock state), light fields entangled to var- 
ious distributions do not interfere with each other due 
to the orthogonality of the Fock states. This reflects the 
which-way information and explains the zero amplitude 
but nonzero photon number in a diffraction minimum. 

If at least one of the modes is a standing wave, the 
angle dependences become much richer. Besides new 
classical maxima given by a± — koxdsmOo ± kixdsinOi, 
the second, "noise," term in Eqs. (|4a|) and (j4b|) is also 
not isotropic. It includes a sum of the geometric coeffi- 
cients squared, which is equivalent to effective doubling 
of the lattice period (or light frequency doubling) lead- 
ing to new features at angles, where classical diffraction 
predicts zero. In Fig. [3l a case of two standing waves 
is shown. Due to the effective period doubling (given by 
2ao,i = 2fco,ia;f^sin0o.i a-nd 2a±), new features at the an- 
gles of, e.g., effective first-order maxima appear, though 
classically only the zero-order maxima are still possible. 

The angle dependence of the photon number variance 
(Ariph)^ determined by (A|Z3p)^ shows anisotropic fea- 
tures due to "period doubling" even for two traveling 
waves. For the coherent state, the light at a maximum 
displays strong noise [(AlDl^)^ = mf^ + 6iV|- + Nk be- 
cause {\D\^) ^ N% + 6A^|. + 7iV| + Nk and = 
N]f + ^k] , stronger than the isotropic component [N]^ 
in highest order of Nk) and new features at 9i = ±7r/2 
(for 6*0 = 0, 2N]^ in highest order of Nk)- In SFm, the 
noise at maxima can be suppressed, while at other angles 
it is nearly equal to that of the coherent state. In MI, 
(A|Dp)^ — 0. Distinguishing between atomic states by 
light statistics is similar to that by the intensity. 

In summary, we have shown that atomic quantum 



states can be nondestructively monitored by measur- 
ing scattered light. In contrast to homogeneous BECs, 
scattering from lattices exhibits advantageous properties: 
suppression of the classical scattering in Bragg minima, 
access to local and nonlocal correlations, angular distri- 
butions richer than classical diffraction. Also, other opti- 
cal phenomena and quantities depending nonlinearly on 
the atom number operators will reflect quantum atom 
statistics [3, 11, [HI, e.g., the dispersion of a medium 



will provide a spectral method of quantum state charac- 
terization |19|]. Exploiting quantum properties of light 
should be applicable to other Bragg spectroscopy setups. 
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